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' ABSTRACT 

7— i ; 

k/ ■ Context. The H.E.S.S. Cherenkov telescope array has been surveying the Galactic Plane for new VHE (>100 GeV) gamma-ray sources. 
rt ■ Aims. We report on a newly detected point-like source, HESS J1943+213. This source coincides with an unidentified hard X-ray source 
^ ' IGR J19443+21 17, which was proposed to have radio and infrared counterparts. 

» I , Methods. We combine new H.E.S.S. , Fermi/LAT and Nancay Radio Telescope observations with pre-existing non-simultaneous multi-wavelength 
. observations of IGR J19443+2117 and discuss the likely source associations as well as the interpretation as an active galactic nucleus, a gamma- 
ray binary or a pulsar wind nebula. 

Results. HESS J1943+213 is detected at the significance level of 7.9cr (post-trials) at RA(J2000) = 19 h 43 ra 55 s ± l' m ± l* ys , DEC(J2000) 
= +21°18'8" + 17^ at ± 20" ys . The source has a soft spectrum with photon index T = 3.1 ± 0.3 stat ± 0.2 sys and a flux above 470 GeV of 
(1.3 ± 0.2 stat ± 0.3 sys ) x 10~ 12 cm -2 s -1 . There is no Fermi/LAT counterpart down to a flux limit of 6 x 10~ 9 cm -2 s -1 in the 0.1-100 GeV 
energy range (95% confidence upper limit calculated for an assumed power-law model with a photon index T = 2.0). The data from radio to VHE 
gamma-rays do not show any significant variability. 

Conclusions. The lack of a massive stellar counteipart disfavors the binary hypothesis, while the soft VHE spectrum would be very unusual in 
case of a pulsar wind nebula. In addition, the distance estimates for Galactic counterparts places them outside of the Milky Way. All available 
observations favor an interpretation as an extreme, high-frequency peaked BL Lac object with a redshift z > 0.14. This would be the first time a 
blazar is detected serendipitously from ground-based VHE observations, and the first VHE AGN detected in the Galactic Plane. 
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1. Introduction 
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The High Energy Stereoscopic System (H.E.S.S. ) is an array of four 
imaging atmospheric Cherenkov telescopes situated in the Khomas 
Highland of Namibia (Aharonian et al. 2006b). H.E.S.S. collaboration 
has been conducting a systematic scan of the Galactic Plane, increasing 
the survey sensitivity and expanding the catalog of very high energy 
(VHE, >100 GeV) sources with new intrinsically faint or distant ob- 
jects. In the survey, which now covers most of the Galactic Plane as 
seen from the Southern Hemis phere, a rich population of VHE galac- 
tic so urces has been discovered ( lAharonian et al.l 2005c. 2006c; Chavesl 
I2009h . 

The majority of known VHE sources are extended beyond the 
~ 0.1° H.E.S.S. point spread function (PSF). The few sources 
in the Galactic Plane that appear point-like are often associated 
with VHE gamma-ray emitting high-mass X-ray binaries (HMXB) 
which include the very well esta blished binaries PSR B 1259-63 
dAharonian et al] l2005bh LS 5 039 dAharonian etaD l2006dh and LS 
I +61 303 jAlbert et al] I2006T) . The point-like VHE source HESS 
J0632+057 is now a strong candid ate for a HMXB system following a 
recen t multi-wavelength campaign (Ahar onian et al.l20 07: Hinto rTet al] 
I2009h . In addition se veral young pulsar wind nebulae (PWN), includ- 
ing the Crab Nebula (Aharonian et al. 2006b), are also unresolved by 
H.E.S.S. . Outside of the Galactic Plane the point-like sources are as- 
sociated with active galactic nuclei (AGN). More than 30 AGN, mostly 
BL Lacertae (BL Lac) objects, have been detected in VHE gamma-rays 
up to nowQ. None have been detected yet serendipitously in the Galactic 
Plane survey. 

Here, we present the discovery of a new unresolved VHE 
gamma-ray source, HESS J1943+213, located in the Galactic Plane 
and spatially coinc ident with an u nidentified INTEGRAL source, 
IGR J19443+2117 dBird et al]|2007T) . This paper is organized as fol- 
lows. In Sect. [2] we present the H.E.S.S. data analysis and results. In 
Sect. [3] we present the sources found within the H.E.S.S. source error 
circle and evaluate the plausibility of each of them being a counterpart 
of the H.E.S.S. source. This section also includes an analysis of archival 
INTEGRAL and Fermi data, as well as of data taken recently by the 
Nancay Radio Telescope. Assuming that all these sources are counter- 
parts to the H.E.S.S. source, Sect. [4] discusses the possible association 
of HESS J 1943+2 13 with an HMXB, PWN and AGN. The conclusions 
are set out in Sect(5] 



2. H.E.S.S. data 

The source was initially discovered in the Galactic Plane scan data col- 
lected between 2005 and 2008. Dedicated observations were performed 
between May and August 2009. These observations were taken in wob- 
ble mode, where the telescopes point in a direction typically at an offset 
of 0.5° from the nominal target position. The offset angle of the scan- 
mode observations varies between 0.5° and 2°. After applying selection 
cuts to the data, to reject periods affected by poor weather conditions 
and hardware problems, the total observation time used for analysis 
amounts to 38.4 hours. This is equivalent to live-time of a 24.8 hours 
on axis due to reduced acceptance for runs (individual continuous ob- 
servation segments of ~ 28 minutes) taken with large offsets from the 
nominal source position. The mean zenith angle of the observations is 
49°. 

The data were calibrated according to lAharonian et al] (120041) . 
Energies are reconstructed takin g the effective optical efficiency into 
account (Ahar onian et al]|2006bl) . Event reconstruction and separation 
of gamma-ray-like events from cosmic-r ay background were made us - 
ing the Model Analysis and standard cuts (de Naurois & Rolland 2009), 
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1 Based on the online catalog for TeV Astronomy TeVCat located at 
http://tevcat.uchicago.edu , provided by Scott Wakely & Deirdre Horan. 
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Fig. 1. Excess map of the field around the position of 
HESS J1943+213, smoothed by the PSF of the instrument (68% 
containment radius of 0.064 deg). 
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Fig. 2. Distribution of squared angular distance (6 2 ) for gamma- 
like events and normalized background events. The dashed line 
denotes a point source profile. 



in which shower images of all triggered telescopes are compared to a 
pre-calculated model by means of a log-likelihood minimization. 

For the spectral analysis, on-source data are taken from a circu- 
lar region of radius 6 2 = 0.01 deg 2 around the source position, and 
the background is subtracted using th e event background rate estimated 
with the reflected background model dAharonian et al.l2006bl) for which 
several off regions at the same distance to the center of the camera as 
the target position are used (excluding the region close to the source). 

A total of No n = 281 on-sources events and A'off = 4086 off-source 
events are measured. The on-off normalization factor is a = ft „/ftoff = 
0.0379, where ft is the solid angle of the respective on- and off-source 
regions. The observed excess is N y = TVon _ aNog = 126 y-rays, cor- 
responding to a significanc e of 8.9o~ pre-tr ials standard deviations ac- 
cording to Equation 6 from lLi & Mai d!983l) . 

In Fig. [T] we plot an excess map of the sky around the 
HESS J1943+213 position smoothed by the PSF of the instrument. The 
angular distribution of events shown in Fig.[2]is consistent with that of a 
point source. After correcting for trials (Aharonian et al. 2006c) the sta- 
tistical significance of the detection is 7.9cr. A fit of a point-like source 
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Fig. 3. Time-averaged VHE spectrum observed from 
HESS J1943+213. The shaded area represents the lcr con- 
fidence level of the fitted spectrum using a power-law hy- 
pothesis. Only the statistical errors are shown and the upper 
limits are 68% confidence level. The lower panel shows the 
residual of the spectrum fit, which are the difference, in each 
reconstructed energy bin, between expected and observed 
number of excess events, normalized to uncertainty on the latter 
((N b s - N exp )/AN exp ). 



convolved with the H.E.S.S. PSF of the excess yields a position at the 
68% confidence level, RA(J2000) = 19 h 43 m 55 s ±l^ t ±l^ ys , DEC(J2000) 
= +21° 18'8" + I7" m ± 20" ys . The source is located in the Galactic plane 
at / = 57.76° and b = -1.29°. The intrinsic source extension, calculated 
at the best fit position is smaller than 2.8' at 3cr confidence level. 

Figure [3] shows the time-averaged differential spectrum which was 
derived using the forward folding technique with an energy thresh- 
old of 470 GeV. Between 470 GeV and ~ 6 TeV this spectrum is 
well described (Log-likelihood chance probability 45%) by a power- 
law dN/dE = a> (£/ lTeVr r with photon index T = 3. 1 ± 0.3 stat ± 0.2 sys 
and normalization at 1 TeV H> (lTeV) = (5.6 ± 0.8 slal ± l.l sys ) x 10~ 13 
cnT 2 s _1 TeV -1 . The 68% confidence level upper limits for the four 
hi ghest-energy bins shown in Fig. [3] were calculated using the method 
of Feldman & Cousins] dl998h . The integrated flux above 470 GeV is 
/(> 470GeV) = (1.3 ± 0.2 stat + 0.3 sys ) x 10~ 12 cnr 2 s"\ which corre- 
spo nds to ~ 2% of the Crab Nebula flux above the same energy thresh- 
old dAharonian et alj|2006bh . In Fig. [4] we show a Lomb-Scargle peri- 
odogram for HESS J1943+213 and in Fig. [5] lightcurves. There is no 
significant variability in the data with chance probabilities of the null 
hypothesis of steady emission respectively of 10.4% and 6.2% for the 
night by night and month by month lightcurve (using statistical uncer- 
tainties only). 

All the H.E.S.S. results presented in this paper have been cross- 
checked with an independent ca libration procedu re and a different back- 
ground discrimination method (lOhm et al.l2009l) . 
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Fig. 4. top: Lomb-Scargle periodogram of HESS J1943+213 
based on the run-by-run lightcurve. The dashed line corresponds 
to the 1% chance probability level. Bottom: distribution of 
power, adjusted by an exponential distribution. In the case of 
a pure white noise, a slope of -1 is expected. 
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Fig. 5. Integral flux /(> ITeV) of HESS J1943+217, measured 
by H.E.S.S. during each observing night (top) and averaged 
per month (bottom). Only the statistical errors are shown. The 
dashed line corresponds in each case to the best fit to a constant 
level. 



3. Counterparts 

We searched for possible H.E.S.S. source counterparts in a number 
of radio, IR, optical and X-ray catalogs. Table [TJ contains the posi- 
tion and spectral information with uncertainties for each source dis- 
cussed in this section by reason of finding them plausible counterpart of 
HESS J1943+213. Fig.[6]shows the H.E.S.S. source position confidence 
contours along with the candidate counterparts positions. 



3. 1 . X-ray counterpart 

The H.E.S.S. source is located within the error circle (4.4') of an 
unidentified hard X-ray INTEGRAL/IBIS source IGR J19443+2117, 
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Table 1. The multi-wavelength fluxes and flux upper limits of potential counterparts of HESS J1943+213. 



Oh^prvatorv 

v_y j v>i vlilui y 


Band 


RA (T20001 


DEC (J2000) 


Position 




a = T - 1 


Flux 


Ref. 






[hh mm ss.s] 


1 




Error 90% 


L J 




[erg cm s ~^ ] 




NVSS (1) 


1 .4 GHz 


19 4^ Sfi 14 


+21 1 


8 26.9 




24.7 


n 57 + n n? '."■>> 


( i 44 + n 0^) y 10" 15 (1 ^> 

^ 1 . ii IT_ \J.\J~J) /\ L\J 


a, b 


NRT (1) 


1 .4 GHz 


19 43 56 


+21 1 


8 26 


15" (8) 


22.8 


0.6 + 0.2 


C1 55 + 0?"l x 10~ 15 (12) 




?MASS (2 > 


9 9 ,, m (K~\ 
z>.z> Liiii yr^-} 


1 4^ Sf, 744 


+21 1 


8 23.38 


< 0.2" 


23.2 




2.37 x 10~ 12 <12a3) 




Swift/UVOT 


0.54 yum (V) 


19 43 56.20 


+21 1 


8 22.95 




22.5 




< 1.6 x 10-' 3 < 12 - 13 > 




ROSAT/HRI (3) 


0.1-2 keV 


19 43 56.2 


+21 1 


8 24 


6" 


23.2 


_(14) 


(3.82 + 0.55) x 10-" <14) 


d 


Chandra^ 


0.3-10 keV 


19 43 56.23 


+21 1 


8 23.6 


0.64" 


23.2 


0.83 + 0.11 


(2.9 + ^)x 10-" 


c 


Swift/XRT <5) 


2-10 keV 


19 43 56.20 


+21 1 


8 22.95 


3.53" 


22.5 


1.04 + 0.12 


(1.83 + 0.04) x 10-" 


f,g 


Swift/BAT (6) 


15-150 keV 


19 43 54.35 


+21 1 


8 08.2 


2.88' 


9.1 




(1.8 + 0.2) x 10-" 


h 


Swift/BAT (5) 


14-195 keV 


19 43 55.2 


+21 1 


8 07.2 




2.9 


i r 3 

1-1 -0.2 


(2.9 + ^)x 10-" 


k 


INTEGRAL/IBIS ,7) 


20-100 keV 


19 44 09.36 


+21 1 


8 25.2 


4.4' 


201 


1.12 + 0.22 


(1.12 + 0.22) x 10-" 


g. i 


INTEGRAL/SPI (7) 


100-200 keV 


19 44 23.13 


+21 16 36.4 




404 




< 2.6 x 10-" 


.i 


Fermi/LAT 


0.1-100 GeV 


19 43 56.23 


+21 18 23.6 




23.2 


1.0 < 15 ' 


< 7 x 10~ 12 




Fermi/LAT 














0.5 (15) 


< 1 x 10~ 12 




H.E.S.S. 


> 470 GeV 


19 43 55 


+21 1 


8 8 


34" 21" 





2.1 ±0.3 slat ±0.2 sys 


a +0.6 sla i+0.5 sys . in-12 





Notes. (1) NVSS J194356+211826; (2) 2MASS J19435624+21 18233; (3) 1RXH J194356.2+211824; (4) CXOU J194356.2+211823; (5) 
SWIFT J1943.5+2120; (6) PBC J1943.9+21 18; (7) IGR J19443+21 17; (8) Error for right ascension only; (9) Angular distance from the centroid 
of H.E.S.S. sou rce position; (10) Spectral index given for all combined radio measurements between 327 and 4850 MHz, listed in SPECFIND V2.0 
dVollmedl2009T) ; (11) Hardness ratio Rate(30-150 keV)/Rate(14-30 keV) = 1.0 ± 0.4; (12) EF{E); (13) The flux is not corrected for extinction; 
(14) Unabsorbed flux obtained from ROSAT/HRI count-rate using a Mission Count Rate Simulator PIMMS, assuming Nu and ax identical with 

Chandra. The flux error includes un certa inties of Ng an d ax ', (15) The spectral ind ex w as assumed for upper limit calc ulation. 

R eferences. (a)ICondon etaild 19981); (b ) JVollmerl d2009l); (c)ISkrutskie et al] d2006l): (d)lROSAT Scientific TeamlfeOOCft: (e ) lTomsick et al.N2009h; 
mlMalizia et all d2007h : (g) lLandi et ai] d2009h : (hi lCusumano et al.l d201Ch : (i) lBird et al.ld2010h : (i) lBouchet et all d2008h : (k) lBaumgartner et all 
d2011h . 



the centroid position of which is about 3.3 arcmin away from the 

H. E.S.S. position. A counterpart to IGR J19443+2117 was detected 
in the soft X-ray d omain by Chandra (CXOU J194356 .2+21 1823; 
Tomsick et all l2009h and Swift (SWIFT J1943.5+2120; lLandi et alj 
2009F with exposure times of 4.8 ks and 11 ks, respectively. The 
same X-ray source was detected in the past by ROSAT/HRI (1R XH 
J194356.2+211824, exposure 1.3 ks. lROSAT Scientific TearnTIOOOl) 

The combined power-law fit to Swift/XRT and INTEGRAL/IBIS 
data gives a spectral index a x = 1.O40 1 and fluxes -F 2 _iokeV = 

I. 83 x 10 ~" ere cnr 2 s" 1 , F 20 -iookev = 1-12 x 10"" erg cm" 2 s" 1 
dLandietal.l |2Q09|). Within the error of the cross-calibration constant 
between XRT and IBIS equal 0.6, the IBIS data is in agreement with 
that of XRT. 

The Chandra observation (performed two years later in 2008), is 
consistent with the Swift results and gives Fo3_| 0ke v = 2.9 x 10~ u 
erg cnr 2 s~' with a spectral index ax = 0.83 (Tomsick et al. 2009]). The 
oldest soft X-ray observation of the source comes from ROSAT/HRI. 
The Mission Count Rate Simulator PIMMS was used in order to obtain 
the flux from the HRI count-rate. Assuming a hydrogen column density 
(see Sect. 13. 2t and spectral index identical to the Chandra source, the 
flux equals F i_2kev = 3.82 x 10~ 11 erg cirr 2 s~', this is in agreement 
with both Chandra and Swift fluxes within errors. 

The INTEGRAL/SPI observations provide an upper limit on the 
source flux above 10 keV of F >I00 keV < 2.6 x 10~" erg cnr 2 s~' 
dBouchet et aD 120081) . An X-ray source consistent with all above 
X-ray observations is also present in the Swift /BAT hard X-ray 
Catalog (PBC J1943.9+21 18. ICusumano et al.ll2010h with Fw-isokeV = 
1.8 x IP'" erg cnr 2 s" 1 , and in the 58 months Swift/BAT catalog 
dBaumgartner et a0l201 lh with Fi4_i95k e v = 2.86 x 10~" erg cirr 2 s~' 
and a x = 1.08. A power -law model fit to the spectrum provided by 
iBaumgartner et all d201 llPI yields ;^ 2 /dof = 2.4/6, whereas a cut-off 



2 a = T — 1; in the following, the subscript of spectral index a de- 
notes a single waveband with R for radio and X for X-ray; or a broad 
waveband with RX for a band extending from radio to X-ray, RO from 
radio to optical and OX from optical to X-ray. 

3 http://heasarc.gsfc.nasa.gov/docs/swift/results/bs58mon/ 



power-law model fit yields ^ 2 /dof = 3/5. There is no evidence for a 
cut-off in the BAT spectrum. 

The public archival INTEGRAL data were analyzed to obtain a 
lightcurve. The IBIS lightcurve was constructed in the 18-60 keV en- 
ergy band, and with individual points corresponding to about ~ 2000 
s (INTEGRAL science window). We used the same methods as thos e 
used in the construction of the 4"' IBIS Catalog dBird et alj 120101) . 
Individual science windows that had exposures of less than 500 s or 
where the source was at a large off-axis angle (> 12° where the in- 
strument flux calibration is not optimized) were excluded from the 
lightcurve. This produced a lightcurve of 1295 observations, spanning 
MJD 52704.16-54604.20, with a cumulative exposure of -2.6 Ms. The 
source is detected at 5.3tr in the 18-60 keV band with a weighted mean 
flux of 0.6 ± 0.1 mCrab; consistent with the reported catalog flux of 
iBird et alj d2010h . The source is weak and its lightcurve is compatible 
with white noise. Fitting the constant to the lightcurve yields ;^ 2 /dof= 
1.2 for 1293 dof. 

Based on the facts that the ROSAT/HRI, Chandra, Swift/XRT and 
INTEGRAL/IBIS fluxes (within the error range of the cross-calibration 
constant) are constant and consistent with each other, that the ROSAT, 
Chandra and Swift sources are only 23" away from the centroid of 
the H.E.S.S. source position, and that there are no other identified X- 
ray sources in the H.E.S.S. error circle, we conclude that these X-ray 
observatories detected the same object and that this object is a very 
likely counterpart to the H.E.S.S. source. Assuming this association is 
correct, in the following, we adopt the Chandra source location with 
its uncertainty for evaluation of potential counterparts in IR and radio 
bands. 

3.2. IR, optical and UV counterpart 

The 2MASS Catalog lists 19 infrared sources within the 90% H.E.S.S. 
error circle, however only one faint unidentified source 2MASS 
J19435624+21 18233 is situated within t he small Chandra er ror circle 
dSkrutskie et"aHl2006l : lLandi et al.ll2009l : lTomsick et al.ll2009h . The an- 
gular separation between the Chandra position and the 2MASS object 
is 0.3". The source is viewed across the Galactic Plane, therefore, its 
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Fig. 6. The dot-dashed lines correspond to (from the innermost) 
68%, 95% and 99% best fit confidence level contours for the 
HESS J1943+213. The positions of the possible counterparts 
are indicated: the partial circle in green/light grey is the (large) 
INTEGRAL error circle for IGR" 19443+2117; the red/grey 
solid circle is radio beam size with NVSS J194356+21 1826 
in center; the blue/grey dashed circle is the NRT position; the 
small magenta/grey solid circle is the ROSAT 90% position 
circle with 1RXH J194356. 2+21 1824 in center; the triangle is 
used for NVSS position; the blue square is used for BAT po- 
sition; the star is used for both CXOU J194356.2+21 1823 and 
2MASS J19435624+21 18233. 



1984: iFurst et al 



3.3. Radio counterpart 

There is only one radio source in the NVSS Catalog within the H.E.S.S. 
error circle (NVSS J194356+211826: ICondon et al.ll998h . It is located 
24.7" away from the best fit H.E.S.S. source position and 3.5" from the 
Chandra source position. Although the Chandra source is outside of 
the NVSS error circle (~ 0.5"), it is still well within the radio PSF 
of 45". The radio source has been detected by various survey pro- 
grams between 3 27 and 4850 MHz a t positions which diffe r from NVSS 
by 17" -40" jTavlor etal.1 119961: Gregory et all 1 1996t iReich et al.l 



--gory e 

19901 : iBecker et al.lll991l : iGregorv & Condon! 1 19911 : 
Griffith et al.ll 199 lb . The sp ectral index ca lculated for nine combined 
measurements is q-r = 0.32 dVol lmer 2009) and t he source flux densit y 
measured by NVSS at 1400 MHz is 102.6 mjy dCondon et al]ll998f) . 
There is no evidence for flux variation over the 12-year time span of the 
collected data. 

Between March and May 2010 the radio source was monitored by 
the Nancay Radio Telescope (NRT), a meridian transit telesc ope with a 
main spherical minor of 300 m x 35 m dTheureau et al .120071) . Thirteen 
observations were performed at two frequencies 1419 and 2392 MHz. 
The data processing consisted of selecting the best drift scans, back- 
ground subtraction, straightening of the baseline by polynomial fit and 
application of a median smoothing. The calibration was corrected by 
the observations of the radio sources 3C123, 3C196, 3C317 and PKS 
0237-23 with the same setup as the main target. 

The average statistical error for all observations is tr slill = 4 mly. 
The systematic errors <r sys = 17 mJy at low frequency and 10 mJy 
at high frequency, take into account the uncertainty in the background 
subtraction and baseline contamination by secondary beams. The for- 
mer mostly affects the low frequency observations due to larger beam 
size (more background sources are found within the beam together with 
the source of interest). The 3-months lightcurve is constant within the 
range of systematic errors. The constant fit to the data yields ^ 2 /dof< 1 
at both frequencies. The average fluxes at low and high frequency are 
111+21 mJy and 86 ± 14 mly, respectively. The average spectral index 
is otr = 0.59 ± 0. 16. The low frequency flux is in agreement with NVSS 
but the spectral index is somewhat higher than the archival observations. 



IR and optical fluxes are strongly affected by extinction. The magni- 
tude lower limits in the J and H bands are 14.53 and 13.36, respec- 
tively, and the source magnitude in the K band is 13.98. The lower limit 
on optical and UV magnitudes of the source from Swift/UVO T obser- 
vations are V > 20.3 and UVM2 > 21.5 iLandi et alj|2009h . Base d 
on maps of Galactic dust infrared emission of ISchlegel eHul fj998l) . 
the extinction coefficients in the direction of the source are A\ = 8.96 
mag, A; = 2.49 mag, A H = 1.59 mag, A K = 1.02 mag for a color ex- 
cess E(B - V) = 2.81 ± 0.25. The extinction corrected magnitudes are 
V > 11.34, J > 12.04, H > 11.77,^0 = 12.96. The exti nction correc- 
tion of the UV magnitude, based on extinction curves of ICardelli et alj 
dl989h . gives U VM2 > 26.04 mag. 

The measured color excess corres ponds to a hydrogen colum n den- 
sity N H = (1.53 ± 0.02) x 10 22 cm" 2 i lPredehl & SclimitIl995T) . which 
is higher than the Nh measure d by the survey of G alactic neutral hy- 
drogen, (8.4 + 1.3) x 10 21 cnT 2 iKalberla et all2005l) . The Nh obtained 
from X-ray Chandra and Swift spect ra are (1.89^ 2 ,) x 10 22 cnT 2 and 
(1.37 +° ! 2 ) x 10 22 cnT 2 respectively llTomsick et al.M2009l: lLandi et al.l 
120091) . thus consistent with the column density obtained from dust 
maps. Therefore, there is no conclusive evidence that the flux of the 
X-ray source is intrinsically absorbed (in spite of observed Nh being 
in excess of Galactic neutral hydrogen measurements). However, dust 
maps used to calcul ate extinction may n ot be precise at Galactic lati- 
tudes lower than 5° (ISchlegel et al.ll 19981) . 

The potential IR counterpart to IGR 119443+2117 was observed 
at ESO by Masetti et al. (private communication, publication in prepa- 
ration), during a campaign aimed at identifying INTEGRAL sources. 
A preliminary IR spectrum of 2 MASS 119435624+211 8233 obtained 
with SofI spectro-imager at NTT (Moorwood et al. 1998) is heavily red- 
dened and does not show any readily apparent emission lines. 



3.4. Fermi/LAT data 

The primary instrument onboard Fermi, the Large Area Telescope 
(LAT) is an electron-positron pair production telescope, featuring solid 
state silicon trackers and caesium iodide calorimeters, sensitive to high- 
energy (HE) gamm a-ray photons from ~ 20 MeV to > 300 GeV 
jAtwood et alJl2009T) . The data were nominally taken in survey mode; 
the observatory is rocked north and south on alternate orbits so that ev- 
ery part of the sky is observed for ~ 30 minutes every 3 hours. 

The LAT analysis dataset spanned MID 54682.7 - 55319.4, cor- 
responding to a period of ~ 21 months. The data were reduced and 
analyzed using the Fermi Science Tools v9rl5 package. The standard 
onboard filt ering, event reconstr uction, and classification were applied 
to the data jAtwood^ t al. 2009), and for this analysis the high-quality 
("diffuse") photon event class is used. Throughout the analysis, the 
"Pass 6 v3 Diffuse" (P6_V3_DIFFUSE) instrument response functions 
are applied. 

Events between 0.1-100 GeV from a 10 degree region of inter- 
est centred on the Chandra X-ray counterpart to HESS 11943+213 
were extracted and a maximum likelihood analysis was performed on 
the data using the gtlike tool. A model was constructed contain- 
ing all the poin t sources within 15 ° that are listed in the Fermi/LAT 
1FGL Catalog dAbdo et alii 20 10a] ); the point sources were modelled 
by a single power-law function and their parameters were frozen to 
those reported in the 1FGL Catalog unless they were within 10° of 
HESS 11943+213. Also incorporated were the currently recommended 
models for the Galactic diffuse emission (gll_iem_vQ2 . fit) and 
isotropic backgrounds; the Galactic diffuse emission is fit with a power- 
law scaling. As HESS 11943+213 does not appear in the 1FGL Catalog 
an additional source was inserted into the model at the Chandra source 
location assuming a power-law spectrum. The likelihood analysis re- 
veals that no significant source is detected at the Chandra source loca- 
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Fig. 7. Spectral energy distribution of the sources discussed in Sectf3] The details of the data and so urces used in this p lot are given 
in Tableland in Sect. I3.2I The dashed black line is the spectral template of elliptical galaxies from lSilva et al.l (Il998l) at a distance 
of 570 Mpc. The dotted (red) line is the blackbody spectrum of a massive star with parameters identical with the companion star in 
the gamma-ray binary LS 5039 at a distance of 25 kpc. The data in IR (2MASS) and visual (Swift/UVOT) bands were corrected for 
extinction. For Fermi/LAT we plot both upper limits for T — 1.5 (dark brown) and T = 2.0 (light brown). 



tion i.e. at the best estimated position of the potential HESS J1943+213 
counterpart. The 95% 0.1-100 GeV flux upper limits calculated using 
a Bayesian approach (Helene 1991) are < 1 x 10" 9 cnr 2 s" 1 for an as- 
sumed spectral index T = 1.5 and < 6 x 10~ 9 cirT 2 s for an assumed 
spectral index T = 2.0. It should be noted that this is a highly complex, 
crowded region in which there is a large contribution from the Galactic 
diffuse emission. 

Both upper limits are shown in Fig. [7] where the limit for T = 1.5 
is plotted in dark brown, and T = 2.0 in light brown. The comparison of 
the Fermi/LAT upper limit with the flux measured by H.E.S.S., implies 
a break in the HE/VHE gamma- ray spectrum located between 100 GeV 
and 500 GeV (or even higher if the spectrum is attenuated by the extra- 
galactic background light), and a hard HE spectrum with T < 2.0. 

4. Discussion 

The positional coincidence between the new VHE gamma-ray source 
HESS J1943+213 and a radio, IR and X-ray source together with the 
lack of any other obvious candidate counteipart strongly suggests an 
association. Figure [7] shows the spectral energy distribution (SED) of 
HESS J1943+213 together with the radio, 2MASS, Chandra, Swift, 
INTEGRAL, and Fermi data. The data are not simultaneous but no sig- 
nificant variability was found within any energy range, therefore, we 
proceed with a contemporaneous discussion of the entire dataset from 
radio to VHE gamma-rays. 

Based on the X-ra y, IR and radio prope rties of IGR J19443+21 17, 
lLandi et all d2009l) and lTomsick et al.l d2009h argued briefly in favour of 
its AGN nature. The detection in VHE gamma-rays helps place tighter 
constraints on the source nature. Because of its point-like appearance 
in VHE gamma-rays, the source is most likely an AGN, a gamma-ray 
binary or a PWN (see Sect. 1). All options are discussed below in light 
of the multiwavelength properties presented in Sect. [3] 

4.1. Gamma-ray binary 

The point-like appearance of HESS J1943+213 in VHE gamma-rays 
and its location in the Galactic Plane may suggest a gamma-ray bi- 
nary. With a handful of known systems, knowledge of these sources as 



a class is still limited. Known gamma-ray binaries are al l radio sources 
with spectral indices a t GHz frequencies «r ~ - 0.5 ijjohnston et al.l 
[1999I ; iRibo et al.|[l999l ; iMassi & Kaufman Bernaddl2009l) . Their X-ray 
spect ra are hard (ax ~ 0-4 - 1.0) with no cut o ffs up to MeV energie s 
(e.g., iTakahashi etafl 120091 ; ISmith et alj|2009l ; lUchivama efaH 120091) . 
The s oftest VHE spectra belong to LS 1+61 303 (T = 2.7 ± 0.5 stat + 
0.2,„ J Albert etal.ll200l) or PSR B1259-63 (T = 2.8 ± 0.2 8tat ± 0.2 sys , 
lAharonian et al.l l 2009t) . In LS 5039 the photon index changes with the 
orbital phase between 1.85 ± 0.25 and 3.09 ± 0.47 ( Aharonian et al. 
I2006al) . The binaries LS 5039, LS 1+61 303 and PSR B 1259-63 have 
been detected by Fermi/LAT and the emission of the former two shows 
a cutoff at a few GeV, suggesting at least two com ponents are require d 
to exp lain their HE and VHE gamma-ray emission dAbdo et al.l2009al lcl 
l2010fh . No HE counterpart has been found for HESS J0632+057. 

The radio, X-ray and VHE spectral properties of HESS J 1943 +2 13 
and its possible counterpart are compatible with what is known to date 
about gamma-ray binaries. The detection of orbital modulation would 
firmly associate HESS J1943+213 with a binary. However the lack 
of modulation in any of the spectral bands does not disqualify the 
HMXB scenario, since for example no modulation has been detected 
in HESS J0632+057. 

All known gamma-ray binaries have a bright, massive star compan- 
ion (O or Be spectral type) dominating the SED in the optical/IR band. 
Identifying the infrared source with a massive star would give credence 
to the identification as a binary. Furthermore, comparing the fluxes to 
those of the massive star in the gamma-ray binary LS 5039 (Fig. |5), 
the IR data are consistent with a massive O star only if its distance is 
> 25 kpc, i.e. if it is located beyond the edge of our Galaxy (radius 
~ 15 kpc) which is unlikely for a massive star. Similarly, comparing to 
the Be counterpart of PSR B 1259-63 (the least luminous and the least 
hot among companions of known gamma-ray binaries) the lower limit 
on the distance is > 22 kpc. A large distance would also imply a 1-10 
keV X-ray luminosity of about 10 36 erg s whereas known gamma-ray 
binary luminosities are about 10 33 - 10 34 erg s" 1 . 

Such a high X-ray luminosity can be explained if the source is 
accreting but, if the source is powered by a pulsar as in PSR B1259- 
63 and based on th e observed correlations in pulsar wind nebulae 
dMattanaetaD 120091) . this would require the spindown power of the 
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pulsar to be higher than that of the Crab: F S 10 38 erg s . Finally, 
the possibility that the massive star is strongly obscured by circumstel- 
lar material, as found in several high-mass X-ray binaries detected by 
INTEGRAL, appears unlikely as the X-ray derived /V H column density 
(see Sect. I3,2I > i s smalle r than what is typically found in these systems 
(e.g. lRahouietaU2008h . 

Despite the general similarity of the radio to X-ray SED between 
HESS J1943+213 and known binaries, the lack of a plausible massive 
stellar counteipart in the optical/IR is the main argument which disfa- 
vors the binary hypothesis. The secondary argument against the binary 
hypothesis is the lack of orbital modulation in any of the bands. 

4.2. Pulsar wind nebula 

When the wind ejected from a rotation-powered pulsar is confined 
by the pressure of the surrounding medium (supernova remnant or 
compressed interstellar gas), a PWN is created. PWNe constitute an 
abundant class of VHE Galactic sources detected by Cherenkov tele- 
scopes . Among them a number are unresolved (e.g ., lAharonian et alj 
l2005al : lDiannati-Atai et al]|2008l : lAcciari et al.ll2010h includ ing one of 
the most powerful, the Crab Nebula ( lAharonian et al.ll2006bl) . 

In general, the SED of PWNe consists of a broad synchrotron 
component extending up to high X-ray energies (~ MeV) and an in- 
verse Compton component at HE and VHE energies (Gaensler & Sland 
120061) . If HESS J1943+213 is a PWN, its gamma-ray (1-30 TeV) to X- 
ray (2-10 keV) flux ratio of 0.04 implies that the nebula is very young 
~ 10 3 y r, and that it contain s a pulsar with a spindown power E ~ 10 38 
erg s -1 (Mattana et al. 2009 ). The high F is also supported by the X-ray 
photon index (Ty ^ 2) dGotthelfll2003l) . 

The estimated PWN age and F is similar to the Crab Nebula. 
Assuming that the luminosity of HESS J1943+213 is comparable to 
the Crab Nebula, we are able to estimate an approximate upper limit 
on the source distance. The low VHE flu x of HESS J19 43+213 being 
only 1.5% of the Crab Nebula (at 2 kpc; lTrimblelll973l) flux, implies 
a distance of < 16 kpc. A comparison with the remaining unresolved 
PWNe, which are also young but less luminous (10-15% of the Crab 
Nebula luminosity), would place HESS J1943+213 at a moderate dis- 
tance of about 6 kpc. 

The young age of HESS J1943+213 in the PWN scenario is 
consistent with it being unresolved for H.E.S.S. . However, all the 
unresolved VHE PWNe are extended in X-ravs dWeisskopf et al.l 
I200d : iNg et alj200"3 : iMatheson & Safi-Harbll2010l : iTemim et al.l2010T) . 
whereas HESS J194 3+213 is point-like as observed by Chandra 
dTomsick et alj|2009l) . Two unidentifi ed H.E.S.S. sources, mo st likely 
associated with bright X-ray pulsars (Aharonian et al. 20063) are em_ 
bedded within very faint diffuse X-ray PWNe, thes e are HESS J1837- 
069 (AX J1838.0-065 5. iGotthelf & HalperrJl2008l) and HESS J1616- 
508 (PSR J1617-5055. lKargaltsev et alj2009f) . A 4.8 ks Chandra obser- 
vation would miss such emission which could explain the point-like ap- 
pearance of HESS 11943+213 in X-rays. However, both of the uniden- 
tified H.E.S.S. sources are extended in VHE gamma-rays and both have 
gamma-ray to X-ray flux ratios greater than unity, much higher than in 
case of HESS J1943+213. 

The hard HE spectrum implied by the upper limits in Fig. [7] (see 
Sect. 13.4b . would be in agreement with the HE spectrum of the Crab 
Nebula , where an IC component has T = 1.64 above 1 GeV (Ab do et alj 
I2010eh . Two other PWN detected by Fermi/LAT (Vela and HESS 
J 1640-465) have soft spectra but these nebulae are older (more evolved) 
than the Crab, and are both extended in VHE gamma- rays. The VHE 
spectrum of HESS 11943+213 is significantly softer than all known 
VHE PWN. The softest PWN have photon indices of 2.7 ± 0.3 and 
there are only a handful of them (Kargaltse v & Pavlovll2010h . This fact 
together with the lack of obvious extended X-ray emission weaken the 
PWN hypothesis. 

4.3. Active galactic nucleus 

Blazars are AGN where a relativistic jet p ointing close to the line 
of sight produces Doppler-boosted emission dBlandford&Reeslll978l : 



Urrv & Padovani 1995). The spectral energy distribution vF v of blazars 
is characterized by two non- thermal components: one peaking in the 
UV to X-ray range, commonly interpreted as synchrotron emission 
from highly relativistic electrons accelerated in the jet, and a second 
peaking in the gamma-ray regime being due to inverse Compton (IC) 
scattering of low-energy photons. Low-energy photons are the syn- 
chrotron radiation emitted either by the same population of electrons 
(synchrotron self-Compton, SSC), or ambient thermal photon (external 
Compton). The Doppler-boosted emission from the jet dominates the 
SED. When blazars are arranged according to their luminosity, the peak 
energy of both components is found to be anti-correlated with luminos- 
ity. High-frequency peak BL Lac (HBL) objects occupy the lowest end 
of the blazar luminosity sequence with the two spectral components 
peaki ng in X-rays (usua lly in the 0.1-1 keV range) and at about 100 
GeV dFossatiet alll998l) . With the IC component peaking at HE, HBL 
objects are the best candidates for VHE sources, and indeed the class 
constitutes the majority of the current population of VHE AGN. 

BL Lac objects are characterized in the optical/IR band by their 
elliptical host galaxy and have featureless optical spectra. Therefore, 
while the available preliminary NTT spectrum does not give a definite 
answer, the apparent lack of emission lines supports a classification of 
HESS 11943+213 as a BL Lac object. 

If the 2MASS IR counterpart of HESS J1943+213 is an elliptical 
host galaxy, its flux ca n be compared w ith a template spectrum of the 
elliptical galaxies, (e.g. lSilva et alj 19981) . The measured flux in K band 
and the flux upper limits in the H, J and V bands, corrected for extinc- 
tion, are consistent with an elliptical galaxy located at a distance > 570 
Mpc (z > 0.14). This is a lower limit on the distance since non-thermal 
emission from the jet in a BL Lac object contributes a large fraction of 
the optical/IR flux. At such distances, absorption of VHE gamma rays 
on the extragalactic background light is expected to significantly soften 
the VHE spectrum. The observed H.E.S.S. spectrum is indeed quite soft 
(photon index of 3.2). The distance also sets the peak X-ray luminosity 
to be > 10 44 erg s" 1 , This is relative ly high for a blazar but still consis- 
tent with spectra of BL Lac objects (Donato et al. 200j]). 

Broad band spectral indices are derived from the observations for 
an assumed host galaxy spectrum and compared to the typical values 
found in blazars. The indices are defined as (tro = - log(Fo/FR)/5.08, 
o-ox = - log(F x /F )/2.6, q-rx = - log(F x /F R )/7.68 where F R , F and 
Fx are the differential fluxes (per unit frequen cy) at 5 GHz, 5000 A, 
and 1 keV respectively dTananbaum et all 19791) . For HESS J1943+213 
the indices are ttRo = 0.36, o-ox = 0.96 and q-rx = 0.56. When com- 
pared to the values of g RO and g ox mea sured for a large group of 
blazars (e.g. Fig. 12 in lPlotkin et atll2010l) . HESS J1943+213 would 
be classified as a radio loud, X-ray strong HBL object. The radio loud 
classification of HESS J1943+213 is also supported by the value of 
\ogR x = log(vL 6cm /L 2 _iokcv) = -3.7, wh ere the radio loud sou r ces ar e 
those with \ogR x > -4.5, as defined by iTerashima & Wilsorj d2003l) . 
The calculated value of log R x is also similar to the average value of 
-3.10 calculated for a large sample of HBL objects, and is significantly 
different than the average -1.27 for low-frequenc y peaked BL Lac ob- 
jects and -0.95 for flat spectrum radio quasars iTerashima & Wilsonl 
120031) . 

ICostamante & Ghisellinil d2002h selected candidate VHE gamma- 
ray emitting BL Lac objects based on their fluxes in radio (F R ) and 
X-ray (Fx) bands measured at 5 GHz and 1 keV, respectively. BL Lac 
objects are expected to be detected at VHE energies if log v R F R > 
-14.8 erg cirT 2 s -1 and logvxFx > -11.4 erg cm -2 s" 1 . This tar- 
get selection criterion has prove d quite successful (e.g. Aharonian et al. 
120081 : lAbramowski et al.ll2010l) . For HESS J1943+213 the fluxes are 
logv R F R = -14.5 and log v x F x = -11.1, which places the source well 
within the group of VHE emitting BL Lac objects. 

The X-ray spectrum of HESS J1943+213 is hard with no ap- 
parent cut-off up to 195 keV. In this respect HESS J1943+213 re- 
sembles extreme BL Lac objects, in which th e synchrotron peak en- 
ergy exceeds 10 keV (Costamante et al. 200l]). The extreme status of 
HESS J1943+213 is also supported by log(F x /F R ) = - 4.3, which is 
highe r than -4.5 as is the case for extreme HBL objects (Rect or et al.l 
120031) . 
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The extreme BL Lac objects are the least luminous blazars in 
the blazar sequence, and yet these are sources in which the dissipa- 
tion of power carried by the jet is most ef ficient and particl e accelera- 
tion mechanisms are close to their limits (Ghisellini 1999). The hard 
X-ray s urveys (e.g. INT EGRAL, Swift) allow selection of extreme 
blazars (Ghisellini 2009). The group of known extreme HBL objects 
includes Mkn 501 1ES 2344+514, PKS 0548-322 and 1ES 1426+428 
dCosta mante et al. 200j|), all of which are VHE sources. The two spec- 
tral components of extreme BL Lac objects peak at energies > 1 keV 
and > 100 GeV for the synchrotron and IC component respectively. The 
HE spectrum in the 100 MeV-100 GeV domain is expected to be hard. 
The lack of an HE counterpart of HESS J1943+213 is unusual. Of all 
the VHE-detected AGN, only fou r (all HBL objects) have not been de- 
tected at GeV energies by Fermi dAb do et al. 201 Obi). The upper limit 
by Fermi is quit e strict but a similar sit uation has been observed e.g., in 
PKS 0548-322 d Aharonian etai]|20ld) . 

Th e Fermi/LPH observations of TeV-selected AGN dAbdo et all 
l2009bl) . show that their HE spectra are hard with photon indices T < 2.0 
where the extreme HBL objects belong to the hardest ones (i.e., 1ES 
2344+514 with T = 1.57 and 1ES 1426+428 T = 1.49). Taking into 
account the FermifLAT upper limit the predicted photon index T < 2 of 
HESS J1943+213 fits well within the HBL hypothesis. 

Blazars are a class of objects known for their variability on all 
timescales. Nevertheless, in the duty cycle of VHE-emitting AGN, there 
are periods without signs of significant flux variability. Accordingly, the 
apparent lack of flux variability in case of HESS J 1943 +2 13 does not 
contradict the hypothesis of it being an HBL object. For instance, HBL 
object PG 1553+113, detected in both HE and VHE gamma rays, has 
shown only weak variability in radio, optical and X-rays and no vari- 
ability in gamma-ra ys despite regular m onitoring of the sky above 200 
MeV by Ferm;/LAT dAbdo et al.|[2010d!) . In addition, recent systematic 
studies of HE variability in b lazars show that HB L objects are the least 
variable subgroup of blazars dAbdo et aTbOlOch . 

Based on the presented arguments, we conclude that the identifica- 
tion of HESS J1943+213 as an extreme HBL object is the most plausi- 
ble hypothesis. 



5. Conclusions 

The H.E.S.S. array of Cherenkov telescopes has detected a new point- 
like VHE gamma-ray source located in the Galactic Plane. The source 
is not variable and has a soft spectrum. The VHE source is coincident 
with the unidentified INTEGRAL source IGR J19443+21 17, which was 
also detected in soft X-rays, radio and IR bands. The positional coinci- 
dence and lack of any other obvious counterpart suggest that these are 
the same object. The available archival data on IGR 119443+21 17 was 
combined with the new H.E.S.S., Fermi and radio observations. This 
allowed us to discuss the source classification as a gamma-ray binary, 
PWN and AGN, the three known types of sources that appear point-like 
in VHE observations. 

The classification as HMXB appears less likely given the absence 
of a massive star. A faint flux of the infrared counterpart attributed to 
a massive star and assuming a similarity of the source SED with the 
known binaries, implies too large a distance (>25 kpc) and in conse- 
quence too high an X-ray luminosity. A secondary less constraining ar- 
gument against the binary hypothesis is the absence of an orbital mod- 
ulation. 

The main argument against the PWN scenario is the soft VHE 
gamma-ray spectrum, which is unusual when compared with the spec- 
tra of a relatively large sample of Te V PWN. In addition, the unresolved 
X-ray emission as observed by Chandra weakens the PWN scenario. 

There are no arguments against the AGN hypothesis. On the con- 
trary, it is strongly supported by the radio, optical and X-ray properties 
being well within the parameter region (set by the broad band spectral 
indices and logarithmic flux ratios) in which high-frequency peaked BL 
Lac objects and TeV AGN are found. The infrared counterpart is com- 
patible with an elliptical galaxy at z > 0.14, which would also explain 
the soft VHE spectrum as attenuation of the intrinsic spectrum on the 
extragalactic background light. The apparent lack of lines in the pre- 



liminary infrared spectrum (if confirmed) could be explained since jet 
continuum emission usually dominates in this energy range. 

HBL objects comprise the majority of AGN detected until now in 
VHE gamma-rays. Statistically, assuming a uniform distribution on the 
sky, one out of 30 AGN can be expected to be situated at Galactic lat- 
itude \b\ < 2°. HESS J1943+213 seems to be one of the most extreme 
HBL objects detected, with a peak frequency well beyond 10 keV and 
a low X-ray to radio ratio. 

The HBL classification would be secured by the identification of an 
elliptical galaxy in infrared images, the measurement of its redshift and 
the detection of correlated variability at radio, optical. X-ray or gamma- 
ray wavelengths. High resolution radio images from VLBI would also 
help to put further constraints on the source identity. 
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